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Abstract 

Minimum-weight triangulation (MWT) is NP-hard. It has a polynomial-time constant- 
factor approximation algorithm, and a variety of effective polynomial-time heuristics that, for 
many instances, can find the exact MWT. Linear programs (LPs) for MWT are well-studied, 
but previously no connection was known between any LP and any approximation algorithm 
or heuristic for MWT. Here we show the first such connections: for an LP formulation due 
to Dantzig et al. (1985): (i) the integrality gap is constant; (ii) given any instance, if the 
aforementioned heuristics find the MWT, then so does the LP. 
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1 Introduction 

In 1979, Garey and Johnson listed minimum-weight triangulation (MWT) as one of a dozen im- 
portant problems not known to be in P nor NP-hard [15]. In 2006 the problem was finally shown 
to be NP-hard [28]. The problem has a sub-exponential time exact algorithm [32], as well as 
a polynomial-time approximation scheme (PTAS) for random inputs |18| . It is still not known 
whether, for some A > 1, finding a A-approximation is NP-hard, but this is unlikely as a quasi- 
polynomial-time approximation scheme exists |31j . MWT has an 0(logn)-approximation algorithm 
|30j . and, most important here, an 0(l)-approximation algorithm called QuasiGreedy [24j. The 
constant in the big-0 upper bound from [2^ is large (we estimate 100,000 or more). 

If restricted to simple polygons, MWT has a well-known 0(n^)-time dynamic-programming 
algorithm jilTl I22j . Polynomial-time algorithms also exist for instances with a constant number of 
"shells" [2] and for instances with only a constant number of vertices in the interior of the region 
R to be triangulated p", §2.5.1], [El 11 |33l [23] . 

Linear program of Dantzig et al. for MWT. Linear programming (LP) methods are one of 
the primary emerging paradigms for the design of approximation algorithms. For many hard combi- 
natorial optimization problems, especially so-called packing and covering problems, the polynomial- 
time approximation algorithm with the best approximation ratio is based on linear programming, 
either via randomized rounding or the primal-dual method. The design of a good approximation 
algorithm is often synonymous with bounding the integrality gap of an underlying LP. 

MWT has several straightforward linear programming (LP) relaxations. Studying their integral- 
ity gaps may lead to better approximation algorithms, or may widen our understanding of general 
methods and their limitations (as standard randomized rounding and primal-dual approaches may 
be insufficient for MWT). 

Dantzig et al. (1985) introduce the following LP (presented here as reformulated by [9j). Below 
A denotes the set of empty trianglesj^ R denotes the region to be triangulated minus the sides of 
triangles in A. The LP asks to assign a non- negative weight Xt to each triangle t G A so that, for 
each point p in the region, the triangles containing it are assigned total weight 1: 

(1.1) minimize c{X) = ^c{t)Xt, subject to X G M>o and (Vp G i?) ^Xt = l. 

teA t3p 

Above, the cost c{t) of triangle t is the sum over the edges e in t of the cost c(e) of the edge, 
defined to be |e|/2 (the length of e), unless e is on the boundary of R, in which case the cost is |e|. 
(Internal edges are discounted by 1/2 since any internal edge occurs in either zero or two triangles 
in any triangulation.) R as specified is infinite, but can easily be restricted to a polynomial-size 
set of points without weakening the LP. (E.g., let R contain, for each possible edge e, two points p 
and g, each on one side of e and very near e.) 

For the simple-polygon case, the above LP finds the exact MWT (every extreme point has 0/1 
coordinates, and so corresponds to a triangulation). This was shown by Dantzig et al. (1985) [3 
Thm. 7], then (apparently independently) by De Loera et al. (1996) [HI Thm. 4.1(i)] and Kirsanov 
(2004) [211 Cor. 3.6.2]. For summaries of these results, see [TUl Ch. 8] and [3S|. Kirsanov describes 

^That is, triangles lying in the region to be triangulated, whose vertices are in the given set of points, but otherwise 
contain none of the given points. 
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an instance (a 13-gon with a point at the center) for which this LP has integrahty gap just above 
1, as well as instances (50 random points equidistant from a center point) that are solved by the 
LP but not by the LMT-skeleton heuristic. 

Other authors have considered edge-based LP^s, mainly for use in branch-and-bound \25 \ [26 l [29 l 
[Ml [3] . These edge-based LPs have unbounded integrality gaps. LPs for maximal independent sets, 
which are well studied, are closely related to all the above LPs, as triangulations can be defined 
as maximal independent sets of triangles (or of edges). The above LPs enforce some, but not all, 
well-studied inequalities for maximal independent sets. 

It is known to be NP-hard to determine whether there exists a triangulation that uses only edges 
in a given subset S [27]. If we change the cost function in the above LP to c{X) = '^Zei^s ^t3e 
the LP will have a zero-cost integer solution iff there is such a triangulation. Unless P=NP, this 
implies that the LP with that cost function has unbounded integrality gapj^ Thus, any bound on 
the integrality gap of the LP with MWT cost function must rely intrinsically on that cost function. 
Similarly, given an arbitrary fractional solution X, it is NP-complete to determine whether there 
is an integer solution in the support of X. These are obstacles to standard randomized-rounding 
methods. 



First new result. We show that LP (1.1) has constant integrality gap. 

This is the first non-trivial upper bound on the integrality gap of any MWT LP. To show 
it, we revisit the analysis of QuasiGreedy (2^, which shows that QuasiGreedy produces a 
triangulation of cost 0(|mwt(G')|), where |mwt(G)| is the length of the MWT of the given instance 
G (and also the cost of the optimal integer solution to the LP). We generalize their arguments to 
show that there exists a triangulation of cost 0{c{X*)), where c{X*) is the cost of the optimal 
fractional solution to the LP. 

Our analysis also reduces the approximation ratio in their analysis by an order of magnitude, 
but the approximation ratio remains a large constant. 



MWT heuristics. Much of the MWT literature concerns polynomial-time heuristics that, given 
an instance, find edges that must be in (or out of) any MWT. Here is a summary. Gilbert observe 
that the shortest potential edge is in every MWT [17| . Yang et al. extend this result by proving that 
an edge xy is in every MWT if, for any edge pq that intersects xy, \xy\ < min{|px|, \py\, \qx\, \qy\} 
P?] . (We refer to the edges satisfying this property as the YXY subgraph.) This subgraph includes 
every edge connecting two mutual nearest neighbors. Keil et al. show (for some /3 > 1) that, if, 
for an edge pq, the two circles of diameter /3 • \pq\ passing through p and q are empty (of other 
vertices), then pq is in every MWT |20j . Cheng et al. strengthen this to /3 ~ 1.17682 [6j. The set 
of such edges is called the fS-skeleton. Das and Joseph show that an edge e cannot be in any MWT 
if both of the two triangles with base e and base angle vr/S contain other vertices [8j . Drysdale et 
al. strengthen this to angle 7r/4.6 [H]. This property of e is called the diamond property. Dickerson 
et al. describe a simple local-minimality property such that, if an edge e lacks the property, the 
edge cannot be in any MWT. Using this, they show that the so-called LMT skeleton must be in 
the MWT dl]. 

^If the LP has bounded integrality gap, it has a zero-cost fractional solution iff it has a zero-cost integer solution. 
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A primary use of the heuristics is to solve some instances of MWT exactly in polynomial time, as 
follows: Given an instance, use the heuristics to identify edges that are in the MWT. If the regions 
left untriangulated by these edges are simple polygons ( equivalently, if the edges connect the given 
points) then find the MWT of each region independently using the standard dynamic programming 
algorithm. (The MWT will be the union of the MWT's of the regions.) According to (1997), 
most random instances with 40,000 points are solvable in this way. 



Second new result. We show that LP (1.1) generalizes these heuristics in that if the heuristics 
solve a given instance as described above, then so does the LP (that is, the extreme points of the 
LP are integer solutions, that is, incidence vectors of optimal triangulations) . 

In fact the LP appears to be stronger than the heuristics, in that some natural instances are 
solved by the LP, but not by the heuristics [21, §3.5] |^ Li this sense, the LP, whose formulation 
requires little explicit geometry, generalizes all of these varied and generally incomparable heuristics. 
This is the first connection we know of between the heuristics and any MWT LP. 

Roughly, the heuristics are based on a combination of (i) local-improvement arguments about 
the MWT and (ii) logical closure (once the heuristic determines the status of one edge with respect 
to the MWT, this in turn determines the status of other edges, and so on). We extend these 
arguments to apply to the optimal fractional triangulation X* . This is possible because (i) X* 
looks "locally" like a MWT and (ii) the LP enforces logical closure of linear constraints on X* . 

After we finished the body of this work, we became aware of and examined additional heuristics 
by Wang et al. [36] and Aichholzer et al. [T]. We conjecture that the LP generalizes them as well. 



An equivalent formulation of the LP. The following constraints are equivalent to the last 



constraints in LP (1.1) (see e.g. [H Thm. l.l(i), Prop. 2.5], [35], or |2H Thm. 3.4.1]) and are useful 



for reasoning about fractional triangulations. For any fractional triangulation X and edge e, 
(1.2) J2 - ^ = [e G boundary(i?)]. 

teleft(e) tGright(e) 

Here left(e) contains the triangles that contain e and lie on one side of e, while right(e) contains the 
triangles that contain e and lie on the other side of e. (If e is on the boundary, take right(e) = 0.) 
The notation [x G S] denotes 1 if x G S and otherwise. 



Practical considerations. Using the O(n^) constraints (1.2) instead of the constraints in (1.1) 
gives an equivalent LP with total size (i.e., non-zeros in the constraint matrix) proportional to the 
number of empty triangles. The empty triangles can be identified, and the LP constructed, in time 
proportional to their number [13]. Their number is always O(n^), but often smaller (e.g. O(n^) in 
expectation for randomly distributed points). 

The time to construct and solve the LP can be further reduced by a preprocessing step based 
on the heuristics — remove any variable Xt if the heuristics prove any edge of t to be out of every 
MWT, and add a constraint J2teMt(e) -^t ~ Steright(e) Xt = 1 if they prove an interior edge e to 
be in every MWT. For randomly distributed points, only 0{n) edges (in expectation) have the 
diamond property, forming O(n^) possible empty triangles, from which the modified LMT skeleton 



^Where G contains the center of a unit circle and n — 1 random points on the circle. 
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can be computed in O(n^) time |1H I12j . On "typical" instances with 10^ — 10^ points, only a 
very small number of variables are left undetermined by the heuristics. (For n random points, the 
expected number is 0(n), but with an apparently astronomically small leading constant [5|.) This 
allows standard LP software to quickly solve the LP, and integer-LP solvers to quickly find the 
MWT. 

Remarks. We do not give an algorithm per se, and the integrality-gap bound, though constant, 
is large. But both results suggest that the LP of Dantzig et al. captures much of the structure of 
MWT. This suggests a clear line of attack for finding an approximation algorithm with reduced 
approximation ratio: study the integrality gap of the LP, trying systematic LP methods such as 
the primal-dual method. If constantly many rounds of lift-and-project (applied to the LP) yield 
an LP with integrality gap 1 + e, this would yield a PTAS. If randomized rounding, primal-dual, 
and similar approaches fail, this will increase our general understanding of the limitations of these 
approaches. 

2 Definitions 

The interior of a segment pq is pq — {p,q}. The interior of a polygon P consists of P minus its 
boundary. Two sets properly intersect (or overlap, or cross) if the intersection of their interiors is 
non-empty. The (Euclidean) length of line segment pq is \pq\- For any set E of segments, \E\ is the 
total length of segments in E. 

A planar straight-line graph (PSLG) is an undirected graph G = {V, E) along with a planar 
embedding that identifies each vertex with a planar point and each edge with the line segment 
connecting its endpoints, so that each edge intersects other edges (and V) only at its endpoints. 
The length of G is the sum of the Euclidean lengths of its edges. G partitions the plane into 
polygonal /acesj^ A face or polygon is empty if its interior contains no vertex. 

A diagonal, or potential edge, of G is any segment pq ^ E connecting two vertices of a face, 
and contained in that face, so that G' = {V, E U {pq}) is still a PSLG. A partition of G is a PSLG 
that extends G by adding (non-crossing) diagonals; equivalently, the faces of the partition refine 
the faces of G. A convex partition of G is a partition whose faces are empty and strictly convex. 
The minimum- length convex partition of G is denoted mcp(G). A triangulation of G is a partition 
whose faces are empty triangles. A fractional triangulation X is a feasible solution to the LP. For 
any potential edge e, the weight of edge e in X, denoted X^, is J2t5e-^t ^ ^he boundary of 
the region to be triangulated, and otherwise half this amount. 

Formally, an instance of MWT is specified by a planar point set V, implicitly defining a PSLG 
G = {V, E) where E contains the edges on the boundary of the convex hull of V. A solution is a 
minimum-length triangulation of G. 

Throughout, we fix an instance G = (y,E) of MWT specified by a given point set V. Unless 
stated otherwise, every graph considered is a partition of G. Since the vertex set V is the same for 
all such graphs, we identify each particular graph by its edge set. 

*Where two points are in the same face if there is a path between them that intersects no edge, with the caveat 
that the term face excludes the single such unbounded region. 
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3 LP (1.1) has constant integrality gap 

Proposition 3.1. Given any instance G = {V,E) of MWT, for any fractional triangulation X, 
there exists an integer solution of value 0{c{X)). Thus, LP (1.1) has constant integrality gap. 

The rest of this section proves the proposition. 

Fix the MWT instance G and an arbitrary fractional triangulation X. For now, also fix an 
arbitrary convex partition CP. (Later, we will specify how to choose CP.) 



Summary of proof. The idea is to define a sort of "rounding" procedure that converts X into 
the desired integer solution. The main step of the procedure converts X into a separate fractional 
triangulation X^ for each face / of CP (covering just /). Next, independently within each face / 
of CP, the procedure replaces the fractional triangulation X^ by the optimal integer triangulation 
of /. The final "rounded" solution is then the union of these integer triangulations (one for each 
face / of CP), of total cost at most X^jgcp ^i-^^) (and, hopefully, 0{c{X))). 

In the second step, since each / is a simple polygon, it follows from known results (e.g. [71 
Thm. 7]; see the introduction) that the cost of the optimal integer triangulation of / is at most the 
cost of X^ . Thus, the integrality gap will be 0(1) as long as the main step triangulates the faces 
so that Y.fC(Xf) = 0{c{X)). 

The proof divides into two parts: (1) defining the rounding procedure and showing that it 
produces a feasible fractional triangulation X^ of each face /, and (2) showing that c{Xf) = 
0{c{X)). 

Proofs of the main theorems are postponed to subsequent sections. 



The rounding procedure. In addition to the given G, X, and CP, let / be an arbitrary face of 
the convex partition CP. 

To convert X into a fractional triangulation X^ of /, start by focusing on just the triangles that 
cross / and have positive weight in X. We "break" each such triangle t into a set t^ of triangles 
within /. Then, in X-^, we give each triangle in t^ weight Xt- 

To break each triangle t into a set t^ of triangles in /, we leverage the concept of edge transposals 
from [24;, (see e.g. Lemma 4.2)]. The reader may skim the details of the definition on first reading. 

Definition 3.1 (transposals of triangles). Given a triangle t crossing face f , the triangulated 
transposal t^ of t (in /) is obtained as follows. First, orienQ each edge e of t so that t lies to the 
right of e. Next, for each edge e of t independently, obtain its edge transposal in f, denoted e^ , as 
follows: 

1. Clip e to e = e n /. 

2. Obtain e^ by sliding each endpoint p of e to an "adjacent" vertex of /: if p is a vertex of 
/, leave it there, otherwise p lies in one edge Y Z of /, slide it to Y or Z, choosing the 
destinations of the endpoints to minimize |e-^| (and breaking ties consistently). 

3. Let e^ inherit e's orientation in the natural way. 
^I.e., order the two endpoints. 
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Now define the (non-triangulated) transposal t-^ of triangle t to be the polygon containing those 
points in f that, for every each edge e oft, lie to the right of its transposal e^ . 

Define the triangulated transposal, t^ , oft to he a minimum-length triangulation of the trans- 
posal t^ . If t^ has no area, then t^ = $. 

The transposal t^ has at most six sides. It might have no area. 

To the right is a collection of trian- 
gles (numbered 1-6, with thin edges), 
blanketing a convex face (with 5 thick 
edges and vertices). To the right of 
that each triangle is clipped into the 
face. To the far right are the (non- 
triangulated) transposals of the triangles. The only triangles whose transposals have positive area 
are triangles 3 and 4; their respective transposals are the 4-gon labeled 3 and the 3-gon labeled 4 
(shown to the far right). (The dashed, gray arrow extending from each edge points towards its edge 




zy. 



transposal. The only edge transposal that does not lie on the boundary of the face / is (xy) 

The fractional triangulation Xf of face / is then obtained (from X) simply by giving weight 
Xt to each triangle t' in the triangulated transposal t^ of i: 

Definition 3.2 (transposal X^ of X in /). The transposal of X in /, denoted X^ , assigns to each 
potential triangle t' in f the weight X^, = 'Yltt'i^tJ -^t- (Here t ranges over triangles that cross f. 
More than one t may contribute to X^,.) 

That Xf is a fractional triangulation of / is not immediate from the definition. This requires 
proof: 

Theorem 3.1. The transposal X^ of X in f is a fractional triangulation of f. That is, it covers 
the points in f uniformly with weight 1. 



Section 3.1 gives the proof. The proof uses two observations — (1) for any given /, the fractional 



triangulation X, restricted to triangles that cross /, can be decomposed into layers, each of which 
looks like an actual triangulation (but possibly extending outside of /), and (2) within each layer, 
the triangles can be collectively "morphed" into their triangulated transposals, while maintaining 
uniform coverage of /. 



Bounding the cost. Having defined the fractional triangulation X^ of each face /, it remains 
to bound their total cost X]/ c{X^) (over all faces of CP). 

The bound will depend on the sensitivity of the edges of CP, defined as follows: 

Definition 3.3 (sensitivity). An edge e is cr-sensitive if, for any potential edge e' that crosses e, 
for each endpoint x of e' , the distance from x to its closest endpoint of e is at most cr\e'\. 

The core of the cost bound is this theorem: 

Theorem 3.2. // every edge in the convex partition CP is a-sensitive, then ^jC^X^) is at most 
3|CP| + 12a c{X). 
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Section 3.2 gives the proof, which is based on several observations. (1) Although a given triangle 
t can cross arbitrarily many faces, and has a transposal t-^ in each of those faces, t crosses at most 
two faces in which its transposal has positive area. Thus, t can contribute to the cost of at most 
two faces. (2) the cost of any transposal t'^ of t (not counting the edges in CP) cannot exceed the 
cost of t by much. This follows from the definition of edge transposals and the sensitivity of CP's 
edges, which imply that, for each edge e of t, the transposal of e cannot be much longer than e. (3) 
Each transposal t^ has at most six sides, so triangulating it to obtain the triangulated transposal 
t-^ increases the cost by a constant factor. 



Thm. 3.2 gives an upper bound of 3|cp| + 12a c{X). 

To use this bound we need CP to have |cp| = 0{c{X)) and a = 0(1). Existing results by 
Levcopoulos and Krznaric get us most of the way there: 

Theorem 3.3 (|24j). For some constant A > 0, and any MWT instance G, there exists a convex 
partition lk of G, whose edges are 4-45'^^''T'S'''t'ive, having total length |lk| < A|mcp(G)|. (Recall 
that MCP(G) is the minimum-length convex partition of G.) 

Proof. Levcopoulos and Krznaric show that what they call the quasi- greedy convex partition has 
these properties: for Property (1), see their Lemma 5.4 and the discussion before it; for Property 
(2), see their Corollary 5.3 |24|. □ 



We now fix the (previously arbitrary) convex partition CP to be the partition LK from Thm. 3.3 
To use the bound in Thm. |3.3[ we need to show that |mcp(G)| is 0{c{X)). 

This is relatively easy. Using the constraints on X and a previous analysis of mcp(G) due to 
Plaisted and Hong [30^ Lemma 10], we show the following bound: 

Lemma 3.1. |mcp(G)| < 18c{X) 

The proof is in Section |3.3[ 

Combining the two theorems and the lemma, the cost of the final integer triangulation is at 
most 

c{Xf) As each / is simple. 

< 3|lk| + 12o-c(X) By Thm. ^ 



< 3A|mcp(G)| +54c(X) By Thm. [S^S 

< 3A- 18c(X) + 54c(X) By Lemma [ST 
= 54(A + 1)c(A:) 



Above A is the constant from Thm. 3.3 (Although A is large, the bound above is still substan- 



tially smaller than the approximation ratio proven for QuasiGreedy in [24J.) 



This (with the proofs of Thm. 3.1, Thm. 3.2, and Lemma |3. 1| below) , proves Proposition 3.1 
3.1 Proof of Thm. ICT 

Let G, X, CP, and / be as above. We start with the observation that the fractional triangulation 
X, restricted to triangles that cross /, can be decomposed into a convex combination of incidence 
vectors of what that we call blankets: 
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Definition 3.4 (blanket). A set B of empty polygons with endpoints in V blankets the face f if 
the union of the polygons contains f and no two of the polygons overlap within f (they may overlap 
outside f). 

If the polygons in B are triangles, the transposal of B (in f), denoted B^ , is the set containing, 
for each triangle t ^ B, the transposal t^ of t. That is, B^ = {t^ \ t G i?}. The triangulated 
transposal, denoted B^ , of B (in f) is just the (multiset) union, over all triangles t ^ B, of the 
triangulated transposal oft. That is, B^ = IJtGB^'^- 

The next lemma says that, over /, X can be decomposed into a convex combination of incidence 
vectors of blankets (containing triangles only). 

Lemma 3.2. There exists a set B of blankets, with a weight eb > for each B & B, such that 
SseB = ^ o-nd, for every triangle t crossing f , Xt = YlseB [t ^ B\eB- 

(Recall "[t G Bf is 1 if t G else 0.) 

Proof. Recall that, for instances consisting of a simple polygon, the LP gives optimal 0/1 solutions 
(e.g., [71 Thm. 7]). We adapt a proof of that property. 

Choose any triangle t that crosses / and has > 0. If any edge e of this triangle crosses (the 
interior of) /, since e has positive weight, there must be a positive- weight triangle t' that has e as 



an edge and lies on e's opposite side (this is implied by Constraint (1.2)). Glue t and t' together to 
form a polygonal region. Continue in this way, growing the polygonal region by repeatedly gluing 
a new triangle to any boundary edge e that crosses /. Stop when the region has no boundary edge 
that crosses /. The triangles glued together in this way must form a blanket B of f. 

Let be the minimum weight of any triangle in B. This gives the first blanket B and its 
weight cb- Subtract from each Xt ioic t £ B . This reduces X's coverage of / uniformly by eb- 
To generate the remaining blankets in B (and their weights), iterate this process as long as X still 
covers / with positive (and necessarily uniform) weight. 

(The process does terminate, as each iteration brings some Xt to zero.) □ 

We will also use the following lemma, whose (long) proof we delay. 
Lemma 3.3. For any blanket B €z B, the triangulated transposal B^ of B triangulates f. 

(The lemma is essentially the theorem we are proving, restricted to the special case when the 
triangles t crossing / have integer weight Xt G {0, 1}, i.e., those with Xt = 1 blanket /.) 

Let /' denote / minus points on potential edges. Fix any point p €z f . We will use the lemmas 
above to show that X^ covers p with weight 1. 

Restrict attention to triangles t that cross /. Recall that X^ is obtained from X by "trans- 
ferring" weight Xt from each triangle t to the triangulated transposal of t. So X^ covers p with 
weight EtEt'ef/b ^ t']^t- 

By Lemma |3.2[ each weight Xt can be split into the sum of the weights of the blankets B 
containing t. That is, Xt = Eseet^ ^ B\eB. 
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Combining these two observations, covers p with weight 

t t'etf B&B 

The final sum on the right, X]('g_g/[P ^ ^']' number of triangles that cover p in the 

triangulated transposal of B. By Lemma |3.3[ this number is 1. Thus, each blanket B contributes 
eB to the coverage of p by X-^ . Thus, X-^ covers p with weight ^^gg e^, which equals 1. 




To finish proving Thm. 3.1, we prove Lemma |3.3[ 

The idea is to morph B continuously into its (non-triangulated) transposal B^ = {t-^ | t G B}. 
Specifically, morph the edges of triangles in B as follows: First, for every triangle edge e, clip e to 
the chord eCif of f, giving a set of chords. Next, for every side YZ of face f (in any order), do the 
following step: simultaneously, for every chord uz having an endpoint z € YZ, slide the endpoint z 
continuously along YZ at unit rate to the corresponding endpoint (Y or Z) ofuz's transposal uz^ . 
As the endpoint z moves, move the chord uz as well (as shown below). 

To the right are the start, middle, and end 
of one step of the morphing process for a single 
side YZ of /. The moving chords are labeled 1- 
9. Chords not touching YZ don't move and aren't 
shown. 

This morphing process clips each edge e of each triangle t & B to the chord e n / , and then morphs 
that chord continuously until it arrives at its transposal, e^^ . We will show below that as the chords 
move no crossings are introduced. Thus, the following invariant is maintained: the regions (each 
one coming from a triangle t £ B) collectively blanket f. 

We will show that each triangle t G B (after being clipped to tCif) is morphed into its transposal 
t-^. Thus, the final set of regions is exactly B-l" , which (by the invariant) must blanket /. 

Since the triangulated transposal B^ of B is obtained from B simply by triangulating each 
polygon in i?-^ (preserving the exact covering of /'), the lemma follows. 

To complete the proof, we show that during morphing no chords cross and that each triangle 
is morphed into its transposal. Consider the step for any side YZ of /. 

Observation: As the endpoint y of a chord ya slides along YZ, the transposal of ya is invariant. 
Proof. By the definition of transposal, (ya)-^ = (za)^ for any z that is (with y) in the interior of 
YZ. Thus, the transposal doesn't change while y stays in the interior of YZ. And, if (e.g.) y is 
moving towards Y, then YU (for some U) is the transposal of yu, so Y must be a closest point in 
{Y, Z} to U. This implies that the transposal of Yu is also YU. Thus, when the chord's endpoint 
y arrives at Y, the transposal of the chord does not change. □ 

The observation implies that the morphing process indeed maps each edge e of each triangle 
t G B to its transposal e-^. 



9 



3.2 Proof of Thm. 



3.2 



3 LP (1.1) HAS CONSTANT INTEGRALITY GAP 




For any two points y, z E YZ, let y ^ z denote that y comes before z if one 
travels from Y afong YZ. Overfoading notation, for any two points u and w on the 
boundary of /, minus YZ, let u ^ w denote that u comes before w if one travels 
from Y to Z along the boundary minus YZ. (In the diagram, y < z and u -< w.) 
Now let yu be an arbitrary chord such that the step slides y towards Y. Let zw be 
an arbitrary chord such that the step slides z towards Z. To finish the proof, we will show that 
y ^ z. Thus, the morphing process does not cause chords to cross. 

Fix U and W such that the transposals of yu and zw are YU and ZW, respectively. 
Observation: Ifw ^ u, then W ^U. (Transposing preserves the order of the non-YZ endpoints.) 
Proof. In the case that u and w both lie in the interior of a single side of /, it must be that the 
transposals of yu and zw are the same (because y and z are also both in the interior of a single 
side, YZ), so U equals W (so W ^ U). In the remaining case (by the definition of transposal), 
there exist two distinct sides e and e' (other than YZ) of / such that u,U £ e and w,W £ e' . Since 



w <u, this implies W <U . 
Observation: U ^ W 

Proof. Since YU is the transposal of yu, point Y must be a closest point in {Y, Z} 
to U; that is, U must lie on the y-side of the bisector of YZ. Likewise W must lie 
on the Z-side of the bisector. It follows from the convexity of / that U ^ W. And, 
since ties are broken consistently in choosing transposals, it cannot be that U = W. 



□ 




□ 



The last two observations imply that u ~< w. Assuming inductively that chords yu and zw are 
non-crossing at the start of the step, this implies that y ^ z. Thus, as y slides towards Y and z 
slides towards Z, the chords remain non-crossing throughout the step. 

□ 



This concludes the proof of Thm. 3.1 



3.2 Proof of Thm. [Ml 

We want to bound the total cost of the fractional triangulations that X induces in all faces / of 
CP, that is, c{Xf). 

In this section, for convenience, we define c{t^ ) = c{tf ) = if t does not cross / or if has 
area zero. We will prove the following lemmas: 

Lemma 13. 4t Any given triangle t crosses at most two faces f in which its transposal t^ has positive 
area. Thus, for a given t, only two faces f have c{t^) > 0. 

Lemma 13. 5t For any triangle t and face f , the cost of t^ minus the edges in CP is at most 2a 
times the cost oft. (Recall that a is the sensitivity of CP's edges.) 

Lemma 13. 6t The cost of the triangulated transposal t-^^ is at most three times the cost of the (non- 
triangulated) transposal t^ . 
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Before we prove the lemmas, we note that they imply the theorem as follows. The total cost is 





E/c(X/) = j:f,tXtc{tf) bydef'nofX/ 

< 3 J2f t -^t ^i^^) by Lemma 3^ 

< 3|CP| + G(TY.tj ^tc{t) [c{tf) > 0] by Lemma [33 

< 3\CP\ + I2a J2t ^tc{t) by Lemma [3^ 
= 3|CP| + l2ac{X). 

Lemma 3.4. Given any triangle t, there are at most two faces f of CP in which t 's transposal t^ 
has positive area. 

Proof. Fix a triangle t = IS.XY Z and consider how the faces of CP can overlap t. Say that a face 
/ is accommodating if t's transposal t^ in / has positive area. 

In the two examples to the right, each dashed 
edge is the edge transposal of an edge of t. Within 
each accommodating face, the (positive-area) trans- 
posal of t is dark. 

We claim that every accommodating face touches 
all three edges of t (where touching an endpoint of 
an edge counts as touching the edge). (For example, 
the accommodating face 2 on the left of the graphic, and 2 and 5 on the right, touch all three edges 
of t. Each other face is non-accommodating and, except for 3 and 4 on the right, touches only two 
edges of t.) 

The claim holds because, if a face / touches only two edges of then ff\t lies within a "corner" 
of t. Then two of t's edges must cross the same two sides (or vertices) of / in the same way, and 
those two edges will have the same transposals (directed oppositely), forcing t^ to have no area. 

Now consider the case that t has a face / that touches the interior of all three edges of t (as in 
the figure to the left, above). Since no other face /' can cross /, no face other than / can touch all 
three edges of t. By the claim, then, only face / might be accommodating, so the lemma holds. 

So assume that no face touches the interiors of all three edges of t. 

By the claim, any accommodating face / still has to touch all three edges of i, but now there is 
at least one edge, say XY , of t whose interior / avoids. Thus, / must touch XY at an endpoint, 
say, Y . (For example, consider the figure on the right above. Faces 2, 3, 4, and 5 touch all three 
edges of t, but not all three interiors.) Since / touches XY at Y , but does not touch the interior of 
Xy, there must be an edge wY of / that extends through the interior of t. Since w is not inside t, 
wY must cut across t to the interior of the edge XZ. Thus, any accommodating face f must share 
some vertex v with t, and an edge of the face must extend from v across the interior of t. 

If there are two accommodating faces, they must extend an edge across t from the same vertex 
V, for otherwise the extending edges would cross inside t. Let this vertex be Y. 

Now consider all edges in CP that extend from Y across the interior of t. Let these edges be 
wiY,W2Y,...,WkY, rotating in order around Y. (In the picture above, k = 3.) CP has k + 1 
corresponding faces fo, fi, ■ ■ ■ , fk, also in order rotating around Y, where and fi share edge 
WiY. By the conclusion of the paragraph before last, only these k+1 faces might be accommodating. 
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To finish, we observe that fi is not accommodating unless i {0,k} (the first or last face). 
Indeed, for i {0, k} edges Wi-iY and WiY of fi extend from Y across t to XZ. Since these edges 
touch at Y, the transposal of XZ in fi is thus just the point Y. Thus, the transposal of t in fi has 



no area. 



□ 





Lemma 3.5. Assume CP's edges are a-sensitive. For any face f and triangle t, the total cost of 
the edges that are in t's transposal in f but not in CP is at most 2a times the cost oft. 

Proof. Let / be any face of CP and e be any edge that crosses /. 

We claim that the length of the edge transposal e^ of e in f is at most 2a times the length of 
e. This claim implies the lemma, because each edge in the transposal of i, but not in CP, is the 
edge transposal e^ of an edge e in t that is not on the 
boundary of /. To finish, we prove the claim. 

For an edge e that crosses a face / one of the 
following cases holds: (1) e is incident to two vertices 
of /, or (2) e is incident to one vertex of / and properly 
intersects one s side of / (as in the left of the figure 
to the right) or (3) e properly intersects two sides s 
and s' of / (as in the right of the figure). 

In case (1), the transposal e^ of e is the same as e, so the claim holds. In case (2), since s is 
o"-sensitive, and e^ is the shortest segment from the endpoint of e to an endpoint of s, \e^ < a\e\. 
In case (3), let Y be an endpoint of e and let Z and Z' respectively be the closest endpoints of 
s and s' to Y . Because e^ is the shortest segment from an endpoint of s to an endpoint of s', 

W\ < \zz'\. 

By the triangle inequality, \ZZ'\ < \YZ\ + \YZ'\. 

Because s and s' are cr-sensitive, \YZ\ and \ YZ'\ are each at most a |e|, proving the lemma. □ 

Lemma 3.6. For any face f and any triangle t, the cost c{t-f) of the triangulated transposal of t 
in f is at most three times the cost c{tf) of the transposal oft in f . 

Proof. The intersection of face / and triangle t is a convex polygon with at most six sides (at most 
three in the triangle boundary, and, alternating with those, at most three from the boundary of 
/). The morphing process described continuously transforms each side of t into its transposal, so 
triangle t is also continuously transformed into its transposal. Thus, t^ is a convex polygon with 
at most six sides. Let set S contain the edges in t^ . Also let set S' contain (up to three) diagonals 
of t^ , connecting alternating vertices around the boundary of t^ such that 5 U 5' partitions into 
a triangulation T (a set of triangles). 

Recall that c(e) = |e|/2 for each e G 5", while c(e) G {|e|/2, |e|} for e £ S. By the choice of 
diagonals, \S'\ < \t^, so c(5') = ^.g^, c(e) < cit^). Clearly c(S) = Eee5c(e) = c(t-^). 

Each edge in S occurs in one triangle in t, while each edge in S' occurs in two. Thus, c(T) = 
EteT<t) = ciS) + 2c{S') <3citf). 

The lemma follows, as c{t^) < c{T). □ 



This concludes the proof of Thm. |3.2 



□ 
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3.3 Proof of Lemma ISTTl 



For every vertex v in the interior of V, define a star at to be a subset of edges incident to v such 
that no two successive edges (around v) are separated by an angle of 180 degrees or more. For 
every vertex v on the boundary of V, define the (only) star at v to consist of the two boundary 
edges incident to v. Let Sram{v) denote the minimum cost of any star at v. Plaisted and Hong 
show |mcp(G)| < 6X;^S'mm(^) |3Ql Lemma 10]. 

We claim 5^in(«) < (3/2) Ee3v Xe\e\. Q 

-^s X]?; '^Ze^v ^e|e| =2 Xe\e\ =2 c{X), this implies the lemma. 

We prove the claim. 

It's easy to see that, for any boundary vertex v, 5'min(w 
just an interior vertex v and its edg es. 

Because X satisfies constraint (1.2), rotating around f, there is a 



X]e9t)^e|e|, so restrict attention to 



sequence ei, 62, . . • , of edges such that each a forms a positive- weight 
triangle with its "neighboring" edge e^+i (ei \i i = k). 

Call this sequence of edges a helix, h. For the rest of the proof, 
let wrap(/i) denote the number of times h wraps around v; let be 
the multiplicity of e in h. By a standard construction the X^s can be 
expressed as a linear combination of incidence vectors of helices. (Similar 
to Lemma 



3.2 




s proof, repeatedly find a helix h, choose weight eh, and 
subtract ehN^ /wT&\>{h) from each Xt, reducing coverage near v by eh.) 
This gives a probability distribution e on helices such that each X^ = Ylh A'^^/wrap(/i). 

Now choose a helix h at random from the probability distribution e. Break (partition) h greedily 
into disjoint groups of contiguous edges such that each group g is maximal subject to the constraint 
that the neighboring edges' angles in g^s total at most 360°. (In the figure, white triangles separate 
groups.) Consideration shows that each group contains a star, and (as neighboring groups are 
separated by at most 180°), there are at least [360wrap(/i)/(360 + 180)] = [2 wrap(/i)/3] groups. 

From the randomly chosen h, choose a group g uniformly at random from /I's first [2wrap(/i)/3] 
groups. For any given edge e, the probability that e is in 5 is at most "^/^eh Nj^ /{2wrap{h)/3) = 
{3/2)Xe. Thus, by linearity of expectation, the expected total length E[\g\] of edges in g is at 
most (3/2) ^gg^ Xe|e|. On the other hand, g contains a star, so E[\g\] > Snim{v). This proves 
Lemma 13.11 □ 



4 LP (1.1) generalizes MWT heuristics 

Fix any MWT instance G = {V,E). It is known that any of the following conditions implies that 
a given potential edge e of G is in every MWT of G. 

/3-skeleton: For f3 w 1.17682, the two disks of diameter f3 \e\ having e as a chord are empty [201 E]. 

yxy-subgraph: Every potential edge pq that crosses e = xy has |e| < min{|px|, \py\, \qx\, \qy\} 
L3iI7j. 

maximality: Every potential edge that crosses e is known to be out 0/ every MWT (see e.g. jllj). 
^Recall that Xe is X^tae ^ "-"^ boundary of the convex huU, and otherwise half this amount. 
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Similarly, any of the following conditions implies that a given potential edge e of G (not on the 
boundary of the region to be triangulated) is out o/ every MWT of G. 

independence: Some potential edge that crosses e is known to be in every MWT e.g. [llj . 
diamond: Neither of the two triangles with base e and base angle it/ 4.6 are empty P 



LMT skeleton: For every two triangles t and t' for which e is locally minimal, one of the edges 
of t or t' is known to be out of every MWT [llj . 

In the LMT-skeleton condition, e is locally minimal for two triangles t and t' ii tCit' = e and t and 
t' together are a minimum- length triangulation of the quadrilateral Q = tL) t' — that is, either Q 
is non-convex, or e is at least as short as the other diagonal of Q. 

Let E* denote the set of edges that can be deduced to be in every MWT by applying the logical 
closure of the above six rules. (Logical closure is necessary because the maximality, independence, 
and LMT-skeleton conditions depend on the known statuses of edges other than e. For example, 
if one of the conditions implies that some edge e is out of every MWT, then the LMT-skeleton 
condition may then in turn imply that some new edge e' is out of every MWT, because e lies on 
one of two triangles t or t' for which e' is locally minimal.) 

Ideally, the set E* gives a partition of G in which every face is empty. If this happens, then 
the remaining edges in the MWT can be found by triangulating each remaining face independently 
using the standard dynamic-programming algorithm, and we say G is solvable via the heuristics. 
According to [llj (1997), most random instances with as many as 40,000 points are solvable via the 
heuristics]^ Next we show that if an instance is solvable via the heuristics, then Linear Program 



(1.1) solves the instance also. 



Proposition 4.1. For any instance G of MWT, let E* be the partition of G defined above. If every 
face of E* is empty, then every optimal extreme point of the LP (for G) is the incidence vector of 
a minimum-length triangulation. 

The remainder of the section gives the proof. The first step is to show that each condition 
above that ensures that an edge is in (or out of) every MWT also ensures that the LP gives the 
edge weight 1 (or 0) in any optimal fractional solution. 

Say that LP ( |1.1[ ) forces a potential edge e to z (where z G {0, 1}) if, for every optimal fractional 
triangulation X* of G, the weight that X* gives to e is z. 

Lemma 4.1. // any of the following conditions holds, the LP forces potential edge e of G to 1. 

1. The f3 -skeleton condition above holds for e. 

2. The Y XY -subgraph condition above holds for e. 

3. (maximality) The LP forces every potential edge that crosses e to 0. 



define the modified LMT-skeleton to be the set of edges that can be deduced to be in every MWT via (the 
logical closure of) just the diamond, LMT-skeleton, maximality, and independence conditions above. The use of 
logical closure is crucial to the effectiveness of the LMT skeleton. 
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Proof idea. Part ([s]) is relatively straightforward: if X* gives weight to every edge that crosses 
e, then no triangle t that crosses e has positive , so the only way X* can cover points near e is 
with triangles that have e as a side. 

The original /3-skeleton and the yxy-subgraph heuristics are shown to be valid for MWT by 
local-improvement arguments: if the condition holds for an edge e that is not in the MWT, then 
a polygon P covering e within the MWT can be retriangulated at lesser cost, contradicting the 
optimality of the MWT [20\ ^ [571 E] • Here we extend those arguments to any optimal fractional 
triangulation X* : if the condition holds and X* does not give e weight 1 , then a polygon P' covering 
e whose triangles have positive weight in X* can be retriangulated (lowering the weight of those 
triangles by e > and raising the weight of other triangles by e), giving a fractional triangulation 
that costs less than X* . 

The original arguments are intricate geometric case analyses, typically taking several pages. 
The arguments do not extend completely to this setting for the following reason: in the MWT 
setting, the polygon P identified for re-triangulation is the union of non-crossing triangles, whereas 
here, in th e fra ctional setting, the polygon P' is the union of triangles that can cross (much as 



in Lemma 



3.2) 





If the triangles in P' don't cross, then the original arguments apply, but in 
general additional analysis is needed. To illustrate, consider the /3-skeleton. Suppose for contra- 
diction that the /3-skeleton condition holds 
for an edge e = uw but it does not occur in 
the MWT. [23 El show that there must be a 
sequence ti,t2, . ■ ■ ,tm of empty triangles in 
the MWT whose union P covers e as shown 
in the left of the figure to the right. Using 
the /3-skeleton condition, they show that this 
union has a triangulation that costs less than 
does ti, . . . , imj contradicting the optimality 
of the MWT. 

In the current context, if e has weight below 1 in X* , then there must (similarly) exist a sequence 
ti,t2, ■ ■ . ,tm of empty triangles with positive weight in X* covering e, but these triangles can cross 
as shown to the right above. We extend their arguments to show that, even if such crossing occurs, 
a triangulation of lower cost can still be found. ooo 

Full proof. Here are the details of the proofs for part[T] (/3-skeleton) and part [2] (yXY-subgraph). 
Part [3] (maximality) is discussed already above, in the proof idea. 

Part [7] ( (3 -skeleton). The original /3-skeleton heuristics are shown to be valid for MWT by local- 
improvement arguments: if an edge e is in the 1.7682-skeleton but not in the MWT, then a polygon 
P covering e within the MWT can be retriangulated at lesser cost, contradicting the optimality of 
the MWT |20| [6] . We briefly sketch their argument and then extend it to any optimal fractional 
triangulation X* . 

Assume for the rest of this section that e goes horizontally from the point a on the left to the 
point b on the right. If ab is not in the MWT, there exists a set of MWT edges that intersect e. 
Let 

ci) • • • ) Cri) be the set of edges indexed in non-decreasing order of their length. If the edges are 
removed from the MWT, an empty polygonal region P results. In [20\ |6] it is shown that P can be 
retriangulated at lesser cost by a set of edges that contains ab. The idea is to generate a sequence 
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of triangulated polygons Pq, . . . ,Pn such that Pq is the degenerate polygon ab, Pn is a triangulation 
of P and Pj-i C Pj. To obtain Pj, Pj-i is expanded to include the endpoints Vj and v'j of ej. 
Assume Vj is above the line through ab and v'j is below it. If both Vj and v'j already lie on the 
boundary of Pj-i then Pj = Pj-i- Otherwise, at least one of them will not be on the boundary 
of Pj-i. Assume without loss of generality Vj is not on the boundary of Pj-i (If v'j is also outside 
Pj-i, it will be dealt with similarly). 

Since Vj is not on the boundary of Pj-i, edge Cj intersects a boundary edge ViVk of Pj-i. 
Consider the sequence 5 of vertices on the path from o to 6 on the boundary of P (there are two 
such paths, but the one above the line through ab is intended). On the sequence 6, vertex Vi is 
the last vertex before Vj that belongs to Pj-i and v^ is the 
first vertex after Vj that belongs to Pj-i- This observation 
allows us to clearly define ViVk in the fractional setting be- 
cause in that setting, polygon P may be self-intersecting and / ' 



may intersect more than one boundary edge of Pj~i in the 




half-space above ab. In general, the triangle AviVjVk contains 
a subsequence 6i of vertices on 6 from Vi to Vj and another ^" 
subsequence 82 from Vj to Vk. The polygon Vi6iVj62Vk is then 
triangulated arbitrarily, and Pj is the union of Pj-i and the 
triangulated polygon Vi6iVj62Vk and possibly another trian- 
gulated polygon to include v'j if v'j is not on the boundary 
of Pj-i- This construction is shown in the right figure. The 
polygon with dashed boundary is P, and Pj^i is a triangu- 
lation of the dark gray polygon. The union of Pj-i and arbitrary triangulations of the light gray 
polygons is Pj which includes the endpoints Vj and v'j of Sj. The light gray polygon above ab is 
Vi6iVj62Vk. The white dots inside triangle AviVjVk are 5i, and the gray dots inside the triangle are 

52- 

|20l |6] show by induction on j that for every j all edges in the triangulation of Pj are shorter 
than ej. By induction all the edges of Pj-i are shorter than Cj-i and \ej-i\ < \ej\. Thus it remains 
to show that all the new edges that are added to Pj-i to form Pj are shorter than Cj. The new 
edges triangulate the polygon Vi6iVj52Vk and they are all contained in triangle AviVjV^, so each 
new edge is shorter than max{|fiT;j|, lujUfcl}. Since ViV^ £ Pj-i, \viVk\ < < It 

remains to show that \viVj\ < \ej\. The argument for VjV^ is similar. The following two facts are 
used for this part of the argument. 

Fact 4.1.1 ([20, Lemma 2]). Let ab be an edge in the 1.17682-skeleton. For any edge pq that 
intersects ab, it holds that \pq\ > max{|a5|, \ ap\, \ aq\, \ bp\, \ bq\}. 

Fact 4.1.2 ([6,j Remote Length Lemma]). Let ab be an edge in the 1.17682-skeleton. Letp,q,r and 
s be four other distinct points of the point set such that pq intersects ab, rs intersects ab, pq does not 
intersect rs, and p and s lie on the same side of the line through ab. Then \qr\ < max{|pg|, \rs\}. 



The argument to show 
maxjlavjl, |a6|} < \ej\ 



< \ej\ is as follows. If Vi lies in triangle Aavjb, then \viVj\ < 



The second inequality holds based on Fact 



4.1.1 



If Vj is outside 



Aavjb, consider the convex hull of the path from a to Vj on Pj. Vertex Vi must lie in some triangle 
AvcVdVj where Vc and Vd are hull vertices. Thus, \viVj\ < max{\vcVj\, {vcVdl, {v^Vjl}. Since Vc and v^ 
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are hull vertices, they were added in the growth process in the past. Thus, the edges Cc and with 
endpoints Vc and Vd were processed before ej , so neither ec nor is longer than ej . Combining this 

< max{|ec|, \ ej\} < \ej\. Using a similar argument, one can 



observation with Fact 



4.1.2 



Jives \VcVj\ 



show that VcVd and v^Vj are both shorter than ej. Thus, \viVj\ < max{\vcVj\, \ vcVd\, \ vdVj\} < \ej\. 
This completes the proof that all the edges of Pj are shorter than ej , and thus the new triangulation 
of P costs less. 

Next we extend the above arguments to any optimal fractional triangulation X* . If X* does 
not give ab weight 1, there is a triangle t with positive X^ that properly intersects ab. For any side 
d of t that intersects ab there must be a triangle t' with positive X^, that has d as a side and lies 



on the other side of d from t. The existence of a triangle t' is a consequence of Constraints (1.2) 
Repeating the same argument for the new triangle(s) gives a set of triangles that cover ab. 

Let r = {ti,t2, . . ■ ,tm) be the sequence 
of triangles in the order they intersect ab in 
the direction from a to b. Triangle ti is in- 
cident on a and is incident on b. All tri- 
angles have a positive weight in X* . The 
triangles in the sequence may or may not 
cross each other as shown in the figure to 
the right. Let P (the shaded area in the fig- 
ure) be the polygon formed by the boundary 
edges of triangles in the sequence. As shown in the right figure polygon P is self-intersecting if 
some of the triangles in the sequence cross. Next, we consider both cases and derive a contradiction 
in each case. 

Case 1. If no triangles in F cross, we can directly apply the technique used in |20) and [6] to 
retriangulate the interior of P at lower cost. Lowering the weight of those triangles in F by e > 
and raising the weight of new triangles by e, gives a fractional triangulation of cost less than c{X*). 





Case 2. If some triangles in F cross, the technique of |20j and [6] cannot be directly applied 
because in the MWT setting, the polygon P identified for retriangulation is the union of non- 
crossing triangles, whereas in this case, P is the union of triangles that cross. 

Let di, . . . ,dn be the set of edges of triangles in F that intersect ab indexed in the order they 
intersect ab (in the direction from a to b). The only part of the argument used in [20] and [6] 
that doesn't go through concerns Fact |4.1^ Fact 4.1.2 holds for any pair of edges pq and rs that 
intersect ab and don't intersect each other. Recall that in the MWT setting the edges intersecting ab 



do not intersect each other, so Fact 4.1.2 holds for any pair of those edges. However, in the current 



case some edges in di, . . . ,dn may intersect each other. Thus, Fact |4.1.2| does not automatically 
hold in this case. 

This issue is resolved by the following technical lemma. 

Sublemma 4.1.1. Let di = pq and dj = rs be two edges of triangles ti,t2, . . . ,tm such that pq 
intersects ab, rs intersects ab, pq intersects rs, and p and s lie on the same side of the line through 
ab. Then ps and qr are both shorter than max{\pq\, \rs\}. 
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Proof. Assume that endpoints p and s are above the Une through ab. Let i be the intersection point 
of di and dj, and assume without loss of generahty that point i is also above the line through ab. 
Let ai and aj be the intersection points of di and dj with ab. Assume i < j so that | aai \ < \ aaj \ . 

We first show that triangle bounded by di, dj and ab 
contains a vertex from the boundary of polygon P. Con- 
sider a maximal subsequence of edges of di, . . . ,dn that 
starts with di and every edge in the subsequence inter- 
sects dj. Let dh = uv be the last edge in this subsequence. 
Note that such dh is well-defined. Let u and v be the end- 
points of dh below and above ab respectively, and let o/j 
be the intersection point of dh and ab. . Clearly ah is 
between Oj and aj. Based on the maximality of the subse- 
quence di, . . . ,dh, edge dh is the last edge that intersects 
dj, so dh+i does not intersect dj. Also dh+i intersects ab 
at a point between ah and aj and shares an endpoint with dh- Edge dh+i cannot be incident on 
V because any edge connecting u to a point on line segment ahaj intersects dj. Therefore, d^+i 
is incident on u. Also, dh+i intersects ahUj, and it doesn't intersect dj or dh- Thus, the other 
endpoint of dh+i must be in the triangle Aaiiaj. Let x be this endpoint. 

The argument so far shows the existence of a point x in 
the triangle bounded by di , dj and ab (see the figure to the 
right). Thus, Zaiiaj < Zaxb < vr/S. The first inequality 
holds because x is in triangle Aaiiaj, and the second in- 
equality holds because ab is in the /3-skeleton (from the def- 
inition of /3-skeleton, it follows that Zawb < vr/S.l for any 
point w in the point set). In triangle Aqir the angle Zqir 
is less than 7r/3, so \qr\ < max{|£(/|, \£r\} < max{|pg|, \rs\}, 
and a similar argument on triangle Apis shows that 
\ps\ < maxi lpql, |rg||. This completes the proof of Sub- 
lemma 



4.L1 




□ 



The rest of the argument is similar to the case 1, where triangles don't cross. The interior 
of polygon P can be retriangulated at lesser cost using the techniques in the original /3-skeleton 
arguments. Finally, lowering the weight of triangles in ti, . . . , f^, by e > and raising the weight 
of new triangles by e, gives a fractional triangulation that costs less than c{X*). 

Part\^ (YXY -subgraph). 

The same argument used for the /3-skeleton holds for the the yxy-subgraph as well. The only 
parts of the argument for the /3-skeleton that uses geometric properties of edges in the /3-skeleton 
are Facts 4.1.1| and 4.1.2 Thus, it suffices to show that the mentioned facts hold for the edges of the 
yxy-subgraph too. Let ab be an edge of the yXF-subgraph and pq be any edge that intersects 
ab. By definition of the yXY-subgraph, \ab\ < min{|pa|, \pb\, \ qa\, \ qb\}, so the union of two disks 
centered at a and b with radius \ab\ doesn't contain p or q. It's easy to show that in triangle Aapq 
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and triangle Abpq the angles Zpaq and Zpbq are both greater than 90°. Hence, 
(4.3) \pq\ > max{|pa|, \pb\, \qa\, \qb\, \ab\}, 



and Fact 4.1.1 holds for the edges of the yxy-subgraph. We next show that Fact 4.1.2 also holds 
for any edge ab of the yXF-subgraph. Let p, q, r and s be four other distinct points of the point 
set such that pq intersects ab, rs intersects ab, pq does not intersect rs, and p and s lie on the 
same side of the line through ab. Let h{q) and h{r) respectively denote the distance of q and r 
from the line through ab. Lemma 2 in |37j states that if h{q) < h{r), then \qr\ < \rs\. Similarly, if 



h{r) < h{q), then \qr\ < \pq\. Hence \qr\ < max{\pq\, \rs\} which proves Fact 4.1.2 for the edges of 



the yxy-subgraph. The rest of the proof is identical to the proof for /3-skeleton edges. □ 

Lemma 4.2. // any of the following conditions holds for a potential edge e of G (not on the 
boundary of the region to be triangulated), the LP forces e to 0. 

1. (independence) The LP forces a potential edge that crosses e to 1. 

2. The diamond condition above holds for e. 

3. (LMT skeleton) For every two triangles t and t' for which e is locally minimal, the LP forces 
one of the edges of t or t' to 0. 

Proof idea. Part ([T]) is straightforward: if potential edges e and e' cross, then the LP covering 
constraint for a point near the intersection of e' and e implies that the total weight of potential 
triangles that have e or e' as sides is at most 1. 

Part the LMT skeleton, is straightforward. If an optimal fractional triangulation X* gives 



e positive weight, then (by constraint (1.2) implied by the LP) there must be two triangles t and 
t' with positive Xj* and X^, whose intersection is e. Edge e must be locally minimal for t and t' 
(otherwise X* could be improved by reducing X* and X^, and raising the weights of the other two 
triangles that triangulate t\Jt'). 

Part ([2]), the diamond condition, is handled as the /3-skeleton and yXF-subgraph are handled 



in the proof idea of Lemma 4.1 



Full proof. As discussed in the proof idea, part [T] (independence) and part |3] (LMT skeleton) are 
straightforward. We give the detailed proof of part [2] (diamond property) here. 

Part^ (diamond property). 

Like /3-skeleton and YXY subgraph, the original diamond heuristic for MWT is proved by local- 
improvement arguments: if the condition holds for an edge e that is in the MWT, then a polygon 
covering e within the MWT can be retriangulated at lesser cost, contradicting the optimality of 
the MWT [El [H] . We first give a summary of the results in [H] and then use them to extend the 
result to any optimal fractional triangulation X* . 
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Suppose e is horizontal and p and q are its endpoints, 
and p is on the left of q. Let Apqw and Apqu be the two 
isosceles triangles with base angle 7r/4.6 above and below 
pq and C be the disk with diameter e as shown in the figure 
to the right. Suppose that Apqw contains a point a' and 
Apqu contains a point b' . If e is in the MWT, a'b' is not 
in the MWT, and there is a set of triangles in the MWT 
that intersect a'b' . Consider the sequence of triangles en- 
countered when tracing a'b' toward a' , starting from edge e 
and stopping with the first triangle that has a vertex inside 
disk C, and let Pi be the polygon formed by the boundary 
edges of triangles in the sequence. Let a be the vertex found 
inside C — if all else fails, then a = a' . In the figure on the previous page. Pi is the dark gray 
area above e. Similarly, consider the sequence of triangles encountered when tracing a'b' toward 6', 
starting from edge e and stopping with the first triangle that has a vertex inside disk C, and let 
P2 be the polygon formed by the boundary edges of triangles in the sequence. In the figure, P2 is 
the dark gray area below e. The boundary edges of Pi are grouped naturally into two chains, one 
from p to a and one from q to o. Vertex a doesn't belong to any of the two chains. Pi is a fan on 
p (or q) if all triangles in Pi are incident on p (or q) . Similarly P2 can be a fan on p (or q) . 

Drysdale et al. in prove the following two facts to show how Pi or P2 or their union can 
be triangulated at lesser cost. 

Fact 4.2.1 (following |14l Lemma 8]). If Pi (or P2) is not a fan, it can be retriangulated at lower 
cost. 

Fact 4.2.2 ([14, Lemma 9]). When both Pi and P2 are fans, then their union can be retriangulated 
at lower cost. 

The above facts contradict the optimality of the MWT. Thus, e cannot be in any MWT. Next 
we extend the above arguments to any optimal fractional triangulation X* . We find polygons P[ 
and P2 corresponding to Pi and P2 in the above argument and show that if X* does not give e 
weight zero, then P{ UP2 can be retriangulated at lesser cost. Lowering the weight of those triangles 
by e > and raising the weight of other triangles by e gives a fractional triangulation that costs 
less than X* . The details of the argument follow. 

If X* does not give pq weight 0, there is a triangle t with positive X^ that has pq as a side. 
Triangle t intersects a'b' . For any side d of t that intersects a'b' there must be a triangle t' with 
positive Xp that has d a side and lies on the other side of d from t. By repeating the same argument 
for the new triangle(s), a sequence F of triangles can be obtained that completely covers a'b'. As 
shown in the figure below triangles in the sequence F may or may not cross each other. The 
left figure shows the case that no triangles cross, while the right figure shows the case that some 
triangles cross. If triangles in the sequence do not cross, the original arguments from [T3] apply. 
However, if triangles in the sequence cross, additional analysis is need. 

®In summarizing the result of [2], we use the same notations and names. The only exception is that [2] uses A 
and B instead of Pi and P2 • 
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Consider the set of triangles in T encountered when tracing a'b' toward a' , starting from edge 
e and stopping with the first triangle that has a vertex inside disk C, and let P[ be the polygon 
formed by the boundary edges of these triangles. Also let a be the vertex found inside C — if all 
else fails, then a = a'. Similarly, define P2 to be the polygon formed by the boundary edges of 
the set of triangles encountered when tracing a'b' from e toward b' until a vertex is inside C. The 
following figure shows P[ and P2 (the dark shaded regions) in two cases. In the left figure there are 
no crossing triangles while in the right figure some triangles in P[ cross. 




we consider the following cases and show how in each case P{ U P2 can be retriangulated at 
lower cost. 

Case 1. P[ is not a fan and some triangles in P[ overlap. 



This case is shown in the figure to the right and 
its magnified version on the next page. In the se- 
quence of triangles in P[ encountered when tracing 
a'b' toward a' , let t' be the first triangle that crosses 
some of the previous triangles in the sequence (the 
light gray triangle labeled in the figure). Let P" be 
the polygon that covers the set of all triangles before 
t' (the dark gray polygon labeled P" in the figure). 
Polygon P" and triangle t' share an edge. Let p'q' be 
this edge such that p' is to left of q' (see the figure on 
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the next page). The boundary edges of P" can be grouped naturally into two chains, one from p 
to p' and one from q io q' . We use the following additional facts from [l^ to show that P" can be 
retriangulated at lower cost. Note that the following facts can be applied to P" because P" is the 
union of triangles that do not cross. 

Fact 4.2.3 ([14, Lemma 7]). // the chain of boundary vertices from p to p' has three consecutive 
vertices x, y and z with \yq\ > \zq\ and the internal angle Zxyz is less than vr, then P" can be 
retriangulated to decrease its cost. (The same is true with q and q' exchanging roles with p and p' , 
respectively.) 

Fact 4.2.4 ([HI Lemma 6]). If the chain of boundary vertices from p to p' has no three consecutive 
vertices x, y, and z with \yq\ > \zq\ that form an internal angle (Zxyz) of less than vr, then 
the clockwise limit on the directions of the boundary edges is perpendicular to qw. (The same is 
true with q and q' exchanging roles with p and p' , respectively, and "counter-clockwise" replacing 
"clockwise" .) 

Let Pi and Pi+i be two consecutive vertices on the chain from p io p' . The clockwise limit on 
the direction of a boundary edge PiPi+i is perpendicular to the line through qw if pj+i is above or 



on the line through pi that is perpendicular to qw. Fact 4.2.3 and Fact 4.2.4 together imply that 



if the clockwise limit on the direction of boundary edges is not perpendicular to the line through 
qw, then P" can be retriangulated to decrease its cost. The rest of proof for Case 1 shows that one 
of the edges on the boundary of P" that t' crosses violates the mentioned limit on the direction of 
boundary edges, and thus, P" can be retriangulated at lower cost. 

In the figure below, the darker shaded area is P" . Triangles in P" don't overlap. Triangle 
t' = Ap'q'v is the first triangle that overlaps some of the previous triangles. Triangle t' and 
polygon P" share edge p'q' , and t' crosses some triangles covered by P" . Thus, t' crosses some 
boundary edges of P" either on the chain from p to p' or on the chain from q to q' . Assume without 
loss of generality that triangle t' crosses an edge on the chain from p to p' . Since edge p'q' is a 
boundary edge of P" , the other two sides of triangle t' (p'v and q'v) should intersect P" . Let piPi+i 
be the first edge on the chain from p to p' that is intersected by q'v when moving from q' to v. 




The line through q'v divides the plane into a half-space above it (the half-space containing point 
w) and a half-space below it. It's easy to see that pi is in the half-space above q'v and Pi+i is below 
it, so pi+i cannot be above the line through pi that is perpendicular to qw. This means that the 
clockwise limit on the direction of the boundary edges is not perpendicular to qw. This combined 
with Fact 4.2.3 and Fact 4.2.4 implies that the interior of P" can be retriangulated at lower cost. 
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Case 2. P[ is not a fan and no two triangles in P[ cross. Since triangles in P[ do not overlap, 
directly implies that P[ can be retriangulated at a lesser cost. 



Fact 



4.2.1 



The previous two cases show that if P[ is not a fan, we can retriangulate some triangles in P[ 
to reduce the cost of triangulation. A similar argument applies to P2, so there remains the case 
where both P[ and P2 are fans. We consider this case next. 

Case 3. P[ and P2 are both fans. In this case, it's easy to see that no two triangle in P[ U P2 can 



cross. Thus, Fact 4.2.2 implies that P[ U P2 can be retriangulated at lower cost. 

In all the above cases, by lowering the weight of some triangles in P{ U by e > and raising 
the weight of some other triangles in the new triangulation by e, a fractional triangulation costing 
less than X* can be obtained, which contradicts the optimality of X* . This completes the proof of 
part [2] (diamond property) and Lemma 4.2 □ 



Assume (as in the statement of Proposition 4.1 ) that the set E* of edges that can be deduced to 



be in every MWT of G gives a partition of G in which every face is empty. It follows from Lemmas 



4.1 and 4.2 (by a simple inductive proof) that every edge that can be deduced to be out of every 
MWT is forced to by the LP, and every edge that can be deduced to be in every MWT is forced to 
1. Thus, in any optimal fractional triangulation X* , no potential triangle t that crosses an edge in 
E* has positive weight X^ . Thus, the optimal fractional triangulations X* are exactly those that, 
for each face / of the partition, induce an optimal fractional triangulation of the simple polygon 
/. It is known (e.g. fTl Thm. 7], [9, Thm. 4.1(i)], [211 Cor. 3.6.2]) that, for any simple polygon /, 
each basic optimal fractional triangulation is the incidence vector of an actual triangulation of /. 
Thus, each optimal extreme point of the LP for G is also the incidence vector of a triangulation of 



G, proving Proposition 4.1 



5 Remarks and open problems 

A main open problem is to find an approximation algorithm with a small constant factor. The bound 
on the integrality gap here, and the bound by Levcopoulos and Krznaric of the approximation ratio 



of their QuasiGreedy algorithm, each include a factor of A (from Thm. 3.3). Their final bound 
includes some additional, comparably large, factors; the analysis here (by focusing on triangles 
instead of edges) simplifies that part of their analysis and avoids those large factors, although it 
does include some (relatively) small constant factors instead. Those small constant factors here can 
clearly be improved somewhat, at the expense of complicating the analysis. But the main challenge 



remaining is to reduce the value of A in Thm. 3.3 



The upper bound shown here on the integrality gap of the LP is constant but quite large. We 
suspect that much better upper bounds can be shown, and that these should lead to an approxi- 
mation algorithm with a better approximation ratio. We suspect that implicit in the analysis here 
is a primal-dual argument; making the dual solution explicit might be a step in this direction. 

The only known lower bounds on the integrality gap are barely above 1. It would be interesting 
to prove a lower bound significantly above 1. 

Since the integrality gap of this LP is above 1, it cannot be used directly to derive a PTAS. 
But applying sufficiently many rounds of lift-and-project to the LP will bring the integrality gap 
to 1 + e. Are only 0(1) rounds required? Does this lead to a PTAS? 
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Docs the LP generalize the heuristics in a stronger sense? Specifically, is there some condition 
(e.g., based on the optimal primal/dual pair) such that, if the condition holds for an edge e, then 
that edge must be in (or out of) every MWT? 

Triangulations optimizing different cost functions (other than the total edge length) are also 
studied in the literature. MWT LPs extend naturally to such problems, by modifying the cost 
function or by restricting the set of empty triangles. (For example, the integrality of the extreme 
points of the LP for the simple-polygon case implies that the simple-polygon result generalizes to 
any linear cost function.) Can results similar to those in this paper be obtained for such problems? 
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